Temperature dependence of optical spectral weights in quarter- filled ladder systems 
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The temperature dependence of the integrated optical conductivity I{T) reflects the changes of 
the kinetic energy as spin and charge correlations develop. It provides a unique way to explore experi- 
mentally the kinetic properties of strongly correlated systems. We calculated I{T) in the frame of a 
t-J-V model for ladder systems, like NaV205, and show that the measured strong T dependence of 
/(T) for NaV205 can be explained by the destruction of short range antiferromagnetic correlations. 
Thus I{T) provides detailed information about super-exchange and magnetic energy scales. 
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For discrete lattice models, which are usually employed 
to describe strong correlation physics in transition metal 
oxides, the integrated optical conductivity (IOC) is di- 
rectly proportional to the kinetic energy [|l|: 



dujaa{uj) (X (-i?kin,a). 



(1) 



which applies for finite temperatures. Eq.(^ is even 
more powerful, since it holds for different polariza- 
tions a individually. This option is particularly use- 
ful for anisotropic systems, like the quarter- filled lad- 
der compounds NaV205 and LiV205 which show highly 
anisotropic optical spectra . Here the T-dependence 
of IOC provides insight into the interplay between kinetic 
energy and interactions. 

Our work is motivated by a recent study of the op- 
tical properties of NaV205 by Presura et al. who 
observed a reduction of I{T) by 12-14% between 4K 
and room temperature. They proposed a fitting formula 
for the optical conductivity (integrated up to 2.25 eV) 
/(T)//(0) = (1 - /exp(-£;o/T)) with / ~ 0.35(0.47) 
and Eq 286(370)K for a{b) polarization, respectively, 
predicting a reduction of almost 50% for 6-polarization at 
several hundred degrees centigrade. As the main absorp- 
tion is near 1 eV one may wonder about the origin of such 
a dramatic change of kinetic energy at low temperatures. 

A key feature of NaV2 05, closely related to its charge 
dynamics ||^ , is the 3D charge ordering transition at 34K 
which is accompanied by the opening of a spin gap . 
The precise structure of the low-T phase is still under dis- 
pute 1^. As discussed below, the charge fiuctuations of 
a single ladder can be mapped onto the Ising model in 
a transverse field (IMTF) ||l^. It has been pointed out 
in Refs. jl^Jll[l that the parameters for NaV205 are such 
that the IMTF is close to its quantum critical point. The 
corresponding soft charge excitations appear at qf, = tt 
and therefore do not directly contribute to <7{uj). Yet it 
was shown in Ref. pd[ | that when including spin degrees 
of freedom the soft charge excitations contribute a small 



absorption continuum in cra{uj) within the charge gap, 
and may explain the anomalous absorption observed by 
Damascelli et al. Presura et al. [|| suggested that 
the low energy excitations at g;, = tt may also explain 
the 30 meV activation energy for charge transport [p^ , 
which is surprisingly low in view of the 0.8 eV optical 
gap. Furthermore it was conjectured js) that these exci- 
tations may cause the T-dependence of IOC's and hence 
Eq in Presura's fit should measure the charge gap. 

The picture which evolves from our calculations is dif- 
ferent. The kinetic energy of the IMTF does not show 
any significant temperature dependence which could be 
attributed to Eq. The temperature scale of the variation 
of the kinetic energy is set by the bare Coulomb inter- 
action. The t-J-V model, however, in agreement with 
the experimental data, displays a dramatic decrease of 
kinetic energy in the range 0.2 J < T < J, where J is the 
exchange integral ]l3| , p^ (order O.leV) of the effective ID 
Heisenberg model which describes the spin dynamics of 
NaV205 g. We attribute the large decrease of I{T) to 
the destruction of short range AF correlations by thermal 
population of local triplet excitations. Thereby IOC's al- 
low to measure the super-exchange contribution to -Bkin- 

At quarter filling, NaV205 can be described by a t- 
J-V model p"5|-p^, as the hopping matrix elements tij 
between vanadium dxy orbitals are small compared with 
the intra-atomic interaction {U ^ i eV): 



H 



YVzjn^rij. (2) 



Elimination of local double occupancies yields the 
Heisenberg term with Jij — itfJU and constrained elec- 



tron creation operators c, 



4.(1- 



which enter 



the densities n,; = '^^r ^1 <j^i,>^ '^^ Coulomb repulsion 
(Vij) between neighbors (Fig.l). Si is the spin-i oper- 
ator at site i. The sums are over spin a =t,i and all 
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neighbor bonds on the trelHs lattice, depicted in 

fig.l. Typical parameters values are included in Fig. 3. 
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FIG. 1. Trellis lattice describing the positions of vanadium 
ions in the (a,b)-plane of Q:'-NaV2 05. 

We begin with the investigation of charge fluctuations 
in the spinless case. Due to the interactions Vij this is 
a nontrivial problem, which can be simplified to a sin- 
gle ladder model without changing the essence of the 
problem, provided Vxy and t^y are not too large. For 
large U/\tij\ and Va/liyl the relevant subspace of one 
electron per rung can be represented by pseudospin op- 
erators Tr, where the eigenvalues ±i of correspond 
to the left/right position of the electron within the nmg 
r, and one obtains the one-dimensional (ID) IMTF pO|: 
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(3) 



Here the Ising interaction, due to Coulomb repulsion, fa- 
vors zigzag charge correlations on the ladder, while the 
kinetic energy i?kin (oc ia) appears as a transverse field 
opposing these correlations. 

The ID IMTF is exactly solvable (l9). The kinetic en- 
ergy, expressed in terms of the dimcnsionless parameter 
h = 2ta/Vb, is given by 



,,cosk + h^ , ,/314e(fc), 

„ , dk — tanhf 

27r /o e(fc ^ 2 



(4) 



where /3 — l/fc^T and e{k) — \/h'^ + 2h cos k + 1 is the 
quasiparticle dispersion which becomes soft at fc = tt as 
h^l (i.e. Eo^VbeiTr)^\2ta-Vb\). 

Fig. 2 shows the temperature dependence of the exact 
results for -Ekin and the pseudo-spin correlation functions 
(CF) Pr = jiYlir'^^r'^^r'+r) ncarcst and next-nearest 
neighbors (r = 1,2). Obviously, the temperature de- 
pendence of i?kin and the CF's is dictated by the bare 
interaction Vb as only T-scale, provided ta ^ Vb . The 
CF's decrease gradually with increasing h = 2ta/Vb, and 
vanish only in the limit h oo. We note one peculiar- 
ity of the exact solution: for h < 1 the magnetization 
(T^) cx £'kin is not strictly monotonic in temperature. 
Here (T^) profits from a decrease of z-correlations. In 
the high temperature limit one finds iS'kin h/T. 



From Fig. 2 we conclude that the charge-only model 
cannot explain the decline of I{T) below room temper- 
ature. To study the problem including spin there are 
two alternatives: (i) the complete spin-pseudospin model 
and (ii) the t-J-V model which we shall choose. 
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FIG. 2. Temperature dependence of the kinetic energy 
(top) and pseudospin CF derived from the exact solution of 
the IMTF for h = 2ta/Vt = 0.5, 1.0, and 1.5. 

For the calculation of the optical spectra and IOC's 
for the t-J-V model we employ, as in Ref. ||ll[], the finite 
temperature Lanczos Technique developed by Jaklic and 
Prelovsek j20|. The exact diagonalization (ED) was per- 
formed for a 4 X 4 site system. As we are studying an 
insulator, we can restrict the calculation of <7ai^) to the 
finite frequency response given by the Kubo formula 

a^iiv) = Re / dre'"^(j„(T)j„), (5) 

^ Jo 

where = i Y.(z,j),a ^y^a (4<tS":<t ~ '^^ ^) 
component of the current operator and i?^ denotes the 
vector connecting sites i and j. In the following, energies 
and temperature are given in eV. 

Figure 3 displays the temperature variation of I{T) 
for the t-J-V model for several sets of parameters on 
a log-log scale. The data show a strong decrease at 
low temperature, which is particularly pronounced for b- 
polarization, while the decrease of the kinetic energy at 
high temperatures (T ^ Vb) is similar to that discussed 
for the pseudo-spin (charge only) model. The low-T vari- 
ation is controlled by the magnetic exchange J oc t1/Vb 
as can be seen by comparing the data for <{, = 0.1 and 0.2 
eV. Furthermore, the results for t^y — and t^y — 0.15 
are similar (for tb = 0.2) ||2l| which is consistent with the 
suppression of super-exchange related to txy^ as discussed 
in Ref. |l|]. 

The near-neighbor spin CF's Sin for tb — 0.2 and 
txy = 0.15 are shown in Fig. 4. The comparison with the 
kinetic energy for these parameters (open circles in Fig. 3) 
clearly illustrates that the low temperature decrease in 
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i^kin is due to the loss of short-range AF spin-correlations, 
whose T- variation is controlled by the magnetic exchange 
J. The largest spin CF's at low temperatures are S14 and 
5i5 consistent with zigzag charge correlations, i.e. large 
Ci4 and C15. As discussed along with Fig. 2, the short- 
range charge correlations, however, exist up to higher 
temperatures, determined by Vb- 
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FIG. 3. IOC for a- (dot-dashed) and 6-polarization (solid 
and dashed lines). Results are given for both t^y = 0.15 
(open) and txy = (filled symbols). Common parameters are 



ta = 0.4, Va^Vb = 0.8 and 



0.9 eV. The set tt = 0.2 



and txy — 0.15 corresponds to spectra shown in Ref. Ill 
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FIG. 4. Spin- Si„ — {SiS^) (top) and charge-correlations 
Cin = (win„) (bottom), (see Fig.l), as function of tempera- 
ture (tt = 0.2, t^y = 0.15). 

Finally, the smallness of C12 below T < 0.4 supports 
the validity of the pseudo-spin model, where only one 
electron per rung is allowed for. We also note that the 
data for Ci„ are consistent with the pseudo-spin correla- 



tion Pr, shown in Fig. 2. 

We emphasize that the kinetic energy is essentially de- 
termined by short-range spin-correlations. In the insu- 
lating state £'fei„ has two major contributions which lead 
to the two distinct T-scales in Fig. 3: (i) valence fluctua- 
tions, which have been discussed already in the context of 
the pseudo-spin model, and (ii) virtual excitations, which 
give rise to magnetic super-exchange whose T-variation 
is controlled by J. 

To understand the low-T behavior of IOC's it is re- 
vealing to study a two-rung system including all terms in 
Eq.(H), except t^y and Vxy, and to employ the spectral 
representation of I{T) (eq.(|l|)): 



E.„ 



En 



(6) 



where the summation is over eigenstates of H and Z is 
the partition function. The eigenstates are characterized 
by {S,pa,Pb), i.e. total spin (S=singlet/triplet) and par- 
ity Pa/b in a/b- direction. The level scheme is sketched 
in Fig. 5. The ground state (s, -f , -I-) is a fully symmetric 
singlet. The lowest excited states are the triplet states 
(t, +, — ) and (t, — , +), where (i, — ) is lower in energy, 
as it is composed of bonding and anti-bonding orbitals 
in a-direction, while (i, — ,-|-) consists of orbitals in b- 
direction and lies already « AOOmeV above the ground 
state. Its thermal population in eq.(^ can be neglected 
at low T . The difference between the two polarizations is 
due to the fact that the current operators ja conserve the 
symmetries {S,Pa,Pb), except for pa —Pa- One finds 
that jb has non vanishing matrix elements between the 
ground state and the first excited singlet with (s, — ), 
whereas the matrix elements involving (<, — ) vanish. 
This has the important consequence that the triplet state 
{t, -f , — ) contributes to Z but not to the numerator of the 
Kubo formula (Eq.(6)), leading to 



Ib{T) « 



AEs 1 -f 36-/^^- 



(7) 



Here the factor 3 reflects the degeneracy of the triplet and 
AEs is the energy difference between the lowest singlet 
states. The situation is different for a-polarization, where 
also (t, — , +) contributes, resulting in an extra term 
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-/3A_Et 



AEs{l + 3e~f>^"t) AEt{3 



(8) 



where k is a ratio of matrix elements. In both cases 
the low-T behavior is governed by local singlet-triplet 
excitations of energy E'^^ and hence by the exchange 
integral J{— E^^) of the effective Id Heisenberg model 
H = J^S^Sr+i ||l^,|l^. The influence is however dif- 
ferent for a- and b-polarization. The T-dependence via 
the optical excitation energies AEs and AEt (~ 1 eV) is 
marginal at low T. 
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We infer that also in extended systems, the low- 
T behavior is due to local singlet-triplet excitations 
with energy, estimated from the two-rung formula, 
E% = Atl/[{Va + Vb)^/TTWJW] for Vb » ta,t,. 
For the standard parameters and % = 0.2, i?°(~70meV. 

As the results of the two-rung system are mainly deter- 
mined by symmetry arguments and energy-scale separa- 
tion, we consider the expressions in Eqs.(7,8) as generic 
results with system-dependent parameters: E^^, k, Ai?s 
and Ai^j. We have determined these parameters by fit- 
ting the ED-data of the 4x4 system (Fig. 3) resulting in 
^ 76 (61) meV for a (b) polarization, respectively, 
and K ^ 0.54. The fit-value for E^^ agrees well with the 
estimate E'^^ . The fit-curves along with the ED-data are 
shown in Fig. 5. 



in origin, and is explained by the destruction of short- 
range spin-correlations as T is increased. We argue there- 
fore, that the observed T-dependence of optical spectral 
weights in NaV205 is explained by the decrease of the 
magnetic super-exchange. This underlines the potential 
of optical experiments to infer the magnetic exchange 
J from the T-dependence of the spectra and to explore 
the changes of the kinetic (and super-exchange) energy 
in correlated systems (particularly in systems with small 
J), which cannot be achieved otherwise. 

It is a pleasure to thank M. Konstantinovic, R. K. Kre- 
mer, Y. Ueda and V. Yushankhai for useful discussions. 
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FIG. 5. IOC for standard parameters with % = 0.2 and 
txy ~ 0.15 eV: fit (solid lines) as described in the text and ED 
from Fig. 3 (symbols). Inset: sketch of optical excitations, of 
the two-rung model. For 6-polarization the triplet transition 
matrix element is zero (dashed). 

A similar fit of the experimental data in yields 
J ~ 39 meV, consistent with values obtained f rom othe r 
experiments, e.g. J ~ 48.2 |^ and 37.9 meV 
According to E'^^ we estimate ti, 



0.15, like 







in [13 


15| 



instead of 0.2 eV used in this study and in Ref. |lll| . 

Finally one may wonder about the change of Ekin at 
the 34 K transition ||^ . As this charge order transition in- 
dicates the onset of 3D long-range order, there is not nec- 
essarily a large change of short range correlations to be 
expected. This view is consistent with the small change 
of /(T) at 34 K in Presura's data. 

In conclusion, we have shown, based on the exact 
CF's of the IMTF, i.e. ignoring spin degrees of free- 
dom, that the low energy excitations associated with 
the (charge) correlations, do not lead to a substantial 
T-dependence of IOC's. However, for the quarter-filled 
t-J-V model we find, in agreement with experimental 
data, a large decrease of IOC's in the temperature range 
0.2 J < T < J. This change in kinetic energy is magnetic 



[2] 
[3] 
[4] 

[5] 

[6] 
[7] 
[8] 
[9] 

[10] 



[11] 

[12] 
[13] 
[14] 



[15] 
[16] 

[17] 
[18] 



[19] 

[20] 
[21] 



[22] 
[23] 



D. Baeriswyl, J. Carmelo, and A. Luther, Phys. Rev. B 
33, 7247 (1986); D. Baeriswyl, C. Gros, and T.M. Rice, 
ibid. B 35, 8391 (1987). 

A. Damascelh et at, Phys. Rev. Lett. 81, 918 (1998). 
S. A. Golubchik et al, J. Phys. Soc. Jpn. 66,4042 (1997). 
M. J. Konstantinovic et ai, Phys. Rev. B 63, 121102(R) 
(2001). 

C. Presura, D. van der Marel, A. Damascelli, and R.K. 

Kremer, Phys. Rev. B 61, 15762 (2000). 

S. Atzkern et al, Phys. Rev. B 63, 165113 (2001). 

T. Ohama et al, Phys. Rev. B 59, 3299 (1998). 

M. Isobe and Y. Ueda, J. Phys. Soc. Jpn. 65, 1178 (1996). 

J. Liidecke et al., Phys. Rev. Lett. 82, 3633 (1999); J.L. 

de Boer et al, Phys. Rev. Lett. 84, 3962 (2000). 

M. V. Mostovoy and D. I. Khomskii, SoUd State Com- 

mun., 113, 159 (2000); M. V. Most ovoy, J. Knoester, and 



:ond-mat/0009464 



D. I. Khomskii, 
M. Cuoco, P. Horsch, and F. Mack, Phys. Rev. B 60, 
R8438 (1999). 

J. Hemberger et ai, Europhys. Lett. 42, 661 (1998). 
P. Horsch and F. Mack, Eur. Phys. J. B 5, 367 (1998). 
For a recent theoretical discussion of super-exchange in 
NaV2 05 see V. Yushankhai and P. Thalmeier, Phys. Rev. 
B 63, 064402 (2001) and references therein. 
H. Smolinski et al, Phys. Rev. Lett. 80, 5164 (1998). 
S. Nishimoto and Y. Ohta, J. Phys. Soc. Jpn. 67, 3679 
(1998). 

J. Riera and D. Poilblanc, Phys. Rev. B 59, 2667 (1999). 
We note that in the quarter-filled case the magnetic prop- 
erties follow from the interplay of H' and 'F' terms in H 
and not from the 'J' term, as Jij = At^j/U is small. 

E. Lieb, T. Schuhz, and D. Mattis, Ann. Phys. 16 407 
(1961); T. H. Niemeijer, Physica 36, 377 (1967); P. 
Pfeuty, Ann. Phys. 57, 79 (1970). 

J. Jaklic and P. Prelovsek, Phys. Rev. B 49, 5065 (1994); 
for a review see Adv. in Physics 49, 1 (2000). 
For small tt = 0.1 this is no longer the case and the 
structure of the ground state is essentially changed and 
determined by t^y 

M. Weiden et al, Z. Phys. B 103, 1 (1997). 

T. Yosihama et al, J. Phys. Soc. Jpn. 67, 744 (1998). 



4 



